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Abstract 

We investigate a class of Cardassian scenarios of the universe evolution in notions of the quali- 
tative theory of dynamical systems. This theory allows us to analyze all solutions for all possible 
initial conditions on the phase plane. In the Cardassian models we find that big-rip singularities 
are present as a typical behavior in the future if n < 0. Some exact solutions for the flat Cardassian 
models as well as a duality relation were found. In turn from the statistical analysis of Knop's SNIa 
data, without any priors on matter content in the model, we obtain that at the 99% confidence 
level this big-rip scenario will reach. The potential function for the Hamiltonian description of 
dynamics is reconstructed from the SNIa data (inverse dynamical problem). We also pointed out 
the statistical analysis results depend oversensitively on the choice of the model parameter M. 
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I. INTRODUCTION 



In this note we apply the qualitative analysis of differential equations to the Cardassian 
models which have become popular during last two years due to some of their interesting 
features. This class of models was proposed as an alternative to the cosmological constant 
m „de, A 3 to explain the enrrent aeeele r atio„ of the Uaive.se Q Q. PVeese and Lewis 

claimed that the Cardassian model, in which the Friedmann-Robertson- Walker (FRW) 
equation is modified, explains the expansion of the universe without any dark energy compo- 
nent. A certain additional term in the FRW equation which may arise from the fundamental 
physics (brane epoch) drives the present acceleration of the Universe. 

The Cardassian model survives several observational tests like the magnitude-redshift for 
the present type la supernovae data z test of the angular size of high- z 

compact radio sources or d^-z test for the SZ/X-ray clusters proposed by Zhu and 
Fujimoto h|. 

The main aim of this paper is to demonstrate the theoretical power in explanation of the 
cosmological problems by the Cardassian model. Moreover, the dynamical system methods 
allow to reveal some unexpected properties of this model. We construct the phase spaces 
for these models and discuss how their structure differs from the canonical model with the 
cosmological constant (Do new types of solutions appear? What is their physical meaning? 
Does any change of parameter values lead to some change of dynamical behavior?). 

From the theoretical point of view it is interesting to analyze all evolutional paths of the 
Cardassian models for all initial conditions. On the other (empirical) hand it is important 
to check which Cardassian evolutional paths fit to astronomical observations, namely type 
la supernovae data. The both issues will be examined in this paper. 

There is a widespread opinion that physically realistic models of the universe should 
possess some kind of structural stability — the existence of too many dramatically different 
mathematical models which agree with observations would be fatal for the empirical method 
of science ^|. A dynamical system is said to be structurally stable if other dynamical sys- 
tems which are close to it (in a metric sense) are topologically equivalent (i.e. modulo 
homeomorphism) . Although the question how to ensure such stability is an open problem 
for the dynamical systems of higher dimension than two j3] i n the case of two-dimensional 
dynamical systems following the Peixoto theorem that states the structurally stable systems 
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form open and dense subsets in the space of all dynamical systems on the plane ^| . More- 
over, in the case of two-dimensional systems there is a simple test of structural stability. 
Namely, if the right-hand sides of the dynamical systems are in the polynomial form, the 
global phase portraits are structurally stable if the number of critical points and limit cycles 
is finite, each point is hyperbolic and there are no trajectories connecting saddle points. 

In the case considered the dynamics is reduced to the form of a two dimensional dynamical 
system with right-hand sides in the polynomial form 

i = ^ = P(x,y), (1) 

y = ^ = Q{x,y), (2) 

where P, Q G C°° class of functions; (x,y) is a differential space Ai, which is useful for 

visualization of the dynamics, called the phase space (or state space). The right-hand sides 

define a vector field V = (P,Q) belonging to the tangent bundle M,. Now we can define a 

phase curve as a integral curve of the vector field. All phase curve with the critical points 

(P(x ,y ) = 0, Q(x ,y ) = 0) constitute the phase portrait of the system. Two phase 

portraits are equivalent if there exists an orientation preserving a homeomorphism which 

maps integral curves of both systems into each other. 

From the physical point of view a critical point represents asymptotic states or equilibria. 

The main aim of qualitative analysis of differential equations is constructing the phase 

portrait of the system. Following the Hartman-Grobman theorem: the nonlinear dynamics 

near the hyperbolic critical points (Vz ReAj ^ 0, where Aj is an eigenvalue of a linearization 

matrix) is qualitatively equivalent to its linear part 

dP dP 
i = —(x ,y ){x - x ) + -Q-{x ,yo)(y - yo) (3) 

y = ^r( x o, yo)(x - x ) + ^-(x 0: yo)(y - yo)- (4) 

Full knowledge of the dynamical system comprises also its behavior at infinity. To achieve 
this one usually transforms the phase plane into a Poincare sphere. Then infinitely distant 
points of the plane are mapped into the equator of the sphere. Of course, the character 
of critical points (which depends on the solutions of characteristic equation A 2 — ATr^l + 
det ^4 = for linearization matrix ^4) is conserved but new critical points can appear at the 
equator. Now, the orthogonal projection of any hemisphere onto the tangent plane gives the 
compactified portrait on compact projective plane. 



There are two main aims of the paper. First, theoretical investigations of the dynamics 
and second, its reconstruction from Knop's SNIa sample (inverse dynamical problem). We 
found that the Cardassian model very well fits the SNIa data and that its unexpected future 
(the big-rip singularity) is consistent with SNIa data on the 99% confidence level without any 
priors on matter content. We have demonstrated how some information about dark energy 
can be deduced from the potential of the Hamiltonian system describing the evolution of the 
universe. Using the potential function instead of the equation of state parameter w = p/p 
seems to be attractive in the context of searching the adequate description of the present 
stage of evolution of the Universe. 

II. QUALITATIVE COSMOLOGY OF THE CARDASSIAN MODEL 

A. Dynamics of Cardassian models filled with single fluid on the (H, p) phase plane 

Let us consider the Friedmann equation generalized by Freese and Lewis [lj to the form 



Let us note that equation (J5J) is the first integral of the generalized Einstein equations 
for the Robertson- Walker metric. The first of these equations is called the Raychaudhuri 
equation and second one is called the conservation equation 



where we assume that the matter content is described by the hydrodynamical energy- 
momentum tensor with energy density p and pressure p. 

For the universe filled with dust (p = 0) system (jUJ) takes the form 



H 2 = § + f(p) 



k 



(5) 



a 



H = -\ (1 + 3/'(p)) (p + p) + i(p + 3/(p)) - H\ 
P = SH(p + p), 



(6) 




P 



3Hp, 



(7) 



where f(p) = Bp n . 
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After the transformation p — > E = ^fp we obtain an equivalent dynamical system 

H = ~E 2 +(l- ^nj BE 2n - H 2 , 

E = -\UE. (8) 

The phase portraits for Cardassian scenarios described by systems (jJJ) and (jHJ) (for different 
values of n) are presented in Fig. ^ 

System (jHJ) has the first integral in the form 

Ik 

PcS = p + Wp n = E 2 + 3BE 2n = 3H 2 + —. (9) 

a 2 

Note that this system is invariant with respect to reflection symmetry E — > —E. As a 
consequence of it we obtain symmetric phase portraits. 

Let us now consider the critical points of system (jHJ at a finite domain of phase space 
{(H,E): E > 0}. They can be of two types: 1) static (H Q = 0, E ^ 0) or 2) non-static 
(H ^ 0, Eq = 0). The static critical point exists if n < 2/3 only. Then 

E = [35(2 - 3n)]^7, # = 0. (10) 

In the opposite case if n > 2/3 there are no static critical points at a finite domain of 
phase space. The non-static critical points must be located on the line {E = 0} (if-axis). 
Therefore in this case n > is required (n > 1/2 for a vector field to be class C 1 ). And we 
can find that it is possible only for H Q = 0. 

From the physical point of view the static critical points represent static universes while 
non-static critical points which are admissible only for n > represent an empty Minkowski 
solution. For discussion of the character of the critical points it is convenient to use the 
linearization matrix at the critical points and the characteristic equation. In the case con- 
sidered we have 

-2H , -\E + n(2- Zr^BE 2 ™- 1 



A 



3 TP 3 fr 

— 2-^0) — 2 



(11) 



Hence for non-static critical points (n > 1/2) we always obtain degenerated critical point 
while for static critical points the characteristic equation takes the form 

A 2 -A(Tr^) + det^ = 0, (12) 
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FIG. 1: The phase portraits for Cardassian scenarios described by the systems of differential 
equations l|7jl.([8|) with different values of parameter n: (a), (b) n = —1; (c), (d) n = 0; (e), 
(f) n = 1. Two types of critical points can exist at a finite domain of the phase space: 1) static 
Hq = and E = Eq only if n < 2/3 or 2) non-static if Hq ^ 0, Eq = 0, hence n = and 
Hq = ±\/S (cosmological constant case). If n > only (0,0) critical point is admissible. The 
static critical points represent the static universe in Cardassian models while stable and unstable 
nodes for n = are de Sitter solutions. The points at infinity can principally be of two types. 
First, if asymptotically standard contribution from the p dominates near the initial singularity, 
and the second one if the Cardassian term dominates the future evolution (future singularity of 
type a — > oo, H — > oo, E — > 0). All these points lie on the intersection of the trajectory of the fiat 
model with a circle at infinity. Future singularities, called big-rip singularities, in the original time 
variable, correspond to a finite value of t = to- We can find it on the phase plane as an intersection 
of a circle at infinity with line {{H, E) : H = ±E n ^fB(2 - 3n)/2, E -> 0}. Note that models with 
n > are structurally unstable because of degenerated critical points at a finite domain of phase 
space. 
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where Tr*4 = and det*4 = XjlE^in — 1). Therefore, if n < 1, then the critical points are 
of a saddle type because the eigenvalues of the linearization matrix are real and opposite 
signs. In the special case n = 1, the static critical point does not exist (see Fig. and non- 
static critical point is degenerated similarly to the case n > 1/2. This means that they are 
corresponding systems because of the presence of non-hyperbolic and degenerated critical 
points in the phase space. 

The behavior of trajectories near the non-degenerated static critical points is equivalent 
to the behavior of linearized system. The linearized system (|SJ) in the vicinity of this point 
is 

H = ±Eo{n-l)[E-Eo], 

E = ~E H, (13) 

where Eq is given by (fTU|) and this point is structurally stable (saddle point) because the 
eigenvalues of the linearization matrix have opposite signs. 
After integration of the system |T3|) we obtain 

H = Cl e~ Xlt + c 2 e X2t , 
E-E = Cl k ie Xlt + c 2 k 2 e X2 \ (14) 

where Ai^ = ±\/ det A = ±E ^ (1 — n)/2 and hi, k 2 are the components of the eigenvectors 
determining main directions along which particular trajectories approach or escape from the 
critical point. 

For completeness of our analysis of the dynamics it is important to investigate the behav- 
ior of trajectories at infinity. For this aim it is useful to compactify M 2 phase space by adding 
a circle at infinity. One can simply do that by introducing of the projective coordinates. 
Two projective maps (z,u), (v,w) cover a circle at infinity given by 

z = 0, -oo < v < oo (H = oo, z = l/H, u = E / H) 

and 

v = 0, — oo < w < oo (E = oo, v = l/E, w — H/E). 
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The system under consideration in the projective coordinates has the form 

d ± = z+ l zu 2 + l -B{3n-2)u 2n z^ 2n \ 
dr o 2 



du 3 

Tr = -2 U + U 



l+ l -u 2 + -B(3n - 2)u 2n z 2 ^ 
6 2 



(15) 



where r is a new time variable such that t — > r: dr = dt/z - it is monotonic function of 
original time variable t only in the domain z > 0. We find that (zq = 0, Mq = 3) (n < 2/3) is 
a critical point of the system |T5|) and for n > 1/2 we obtain additionally the critical point 
(zq — 0, Uq — 0). In this case we have 



Tr ^ = 2 + 3 W °' d6t ^ = 12 ( 6 + M °) ( M ° ~ 1 



2 3 u 12 V 

After simple calculation we obtain that the critical points (0, ±v3) are stable nodes located 
on a circle at infinity as intersection points with line Eq = ±Hq. The eigenvalues of the 
linearization matrix at the point (zq, u = ±v3) are Ai = 1 and A2 = 3/2 i.e., they are real 
and positive therefore the critical points are unstable (TrA = 5/2 > 0) for both Hq > 
and Hq < if considered in the original time parameterization (see Fig. [IJ. The dynamical 
system methods offer the possibility of a simple detection of some other unexpected critical 
points in ACDM model. Note that the critical points at infinity can also appear on the 
intersection of a circle at infinity with lines {(H,E): H = ±E n y^B(2 — 3n)/2}. Then 
E — ► and H — ► 00 only if n assumes negative values. This critical point represents an 
asymptotic state of the evolution dominated by the Cardassian term as E — > (a — > 00). 
As a result of such domination for n < 2/3 we obtain dangerous big- rip singularities. From 
Fig. 1. one can conclude that if n < they are typical asymptotic stages of the evolution 
in the future. In this stage H is going to infinity while energy E is reaching zero. This 
critical point corresponds to r = In a = 00 or a finite value of the cosmological time t. In 
the projective coordinates at this point we have 



^^ ± fEf, o = . (16) 

The phase plane (H, p) can also be useful in analysis of model under consideration in 
which the effects of bulk viscosity are included. They are equivalent to introducing the 
effective pressure p c s = p — 3£H with exact dependence on the Hubble function; here £ is 
called viscosity coefficient and in the Belinskii parameterization takes the form £ = ap m . 
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For this case we can rewrite system (JBJ) in the form 

Bp n - H 



3 -n(l-3ap m - 1 H) 



2 

p = -3H(p-3ap m H). (17) 

The phase portraits for representative cases are demonstrated in Fig. |21 

As it is well known there are some major motivations to study such dark energy models 
which have de Sitter solution as a global attractor at late times Q| . From Fig. |2 we can 
observe the existence of de Sitter solution as a global attractor. Note that there are both 
open and closed universes approaching this critical point. All such critical points are the 
solutions of the following system 

p cff - 3H 2 = or p + 3Bp n = 3H 2 , 

p = 3£(p)H = 3ap m H. (18) 

The analysis of the evolution of FRW cosmological models with dissipation, by using (H, p) 
plane, can be performed for the general form of quintessence coefficient w = w(p, H), p = wp, 
or for the models with causal viscosity }lfil |. 

While the evolution of Cardassian model of the universe filled with single fluid, for which 
the equation of state assumes the form p = Wip, is described by the system of differential 
equations in the phase variables (H,p), it is useful to consider a different representation of 
the dynamics basing on the Hamilton methods. 



B. Particle like description of Cardassian models 

In our previous works it was shown that the dynamics of any quintessential FRW 

cosmologies can be reduced to the form of a one-dimensional Hamiltonian flow. In this 
approach we postulate that the equation of state has the general form p = w(a(z))p in 
which the equation of state factor w can be parameterized by the scale factor a or redshift 
z. We obtain that any FRW dynamics with such a fluid is described by the system of 
differential equations 

x = y, 
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FIG. 2: The phase portraits for Cardassian scenarios described by the system of differential equa- 
tions l)17|) for B = a = 1 and for different values of parameters (n, m). There are critical points of 
two types: static which lie on the p-axis and non-static which are the intersection of the trajectory 
of the flat model (dashed line) with the boundary of the weak energy condition p = 3ap m H (dotted 
line). Note that obtained phase portraits are structurally stable at a finite domain because of the 
presence of nodes and saddles. 



where x = a/ do, a is the scale factor of the universe, and the Hamiltonian of the system 
under consideration assumes 



n = j + v(x) = o. 

The potential function V(x) for system ()19j) takes the form 

,.2 



V(x) 



6 



(20) 



(21) 
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For standard Cardassian model of the universe, filled with matter, for which the equation 
of state is p = Wip (wi = const.), we can simply obtain that 



V(x) 



1 



-1 — 3u>; 



^Card,0^ 



-1-37 



(22) 



where 7 = n(l + u?j) — 1. 

The advantage of representing dynamics in terms of Hamiltonian is the possibility to 
discuss the stability of critical points only from the geometry of the potential function. In 
the general case of the system with natural Lagrangian function, if a diagram of the potential 
function has a maximum and is upper convex, then it corresponds to the stable attractor. 
In our case the Hamiltonian function takes the special form of a simple mechanical system. 
Then the linearization matrix A is such that: 

'd 2 V 



TiA = 0, and det A 



dx 1 



(23) 



{x=xo,y=Q) 



Therefore from the characteristic equation A 2 — A(Tr^4)+det A = we obtain that eigenvalues 
of the linearization matrix should be 1) real and opposite signs or 2) purely imaginary and 
conjugated. Therefore the possible critical points in a finite domain of the phase space are 
only saddles (for the first case) or centers (for the second case). In the case considered we 
have 



1 + 37 Dcard 
1 + 3Wi Qi n 



3(7-™;) 



det A 



|x - 3(1+ ^l\ (l + Wi)(n + !)(! + 3^). 



(24) 



Therefore u>; = — 1, —1/3; n = — 1 are always bifurcations values of model parameters for 
which the global structure of the phase space changes qualitatively. If critical point exists 
((1 + 37) (1 + 3wi) < 0) then on the phase plane (x, x) we obtain saddles at a finite domain 
if only (1 + Wi)(l + 3wi)(n + 1) > 0. In the opposite case we obtain centers which are 
structurally unstable. 

Now, let us concentrate on the models for which the potential function is upper convex. 

nnn 

As we demonstrate in our previous papers ( [7|, [17|, |18[ ) such a form is favored by SNIa data 
(see also next section with reconstruction of V(a) from Knop's sample of SNIa data). It is 
interesting that classification of all possible (representative) evolutional paths in the phase 
plane can be simply performed in terms of a pair (w,n) (see Fig|3J). There are presented 
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domains (shaded regions) on the plane (w, n) for which the Cardassian models, described 
by the potential function (|22j). start to accelerate for a > in Fig. El In these cases the 
potential function V(x) of the system has a maximum for a > and is upper convex. The 
phase portraits for the representative cases of Cardassian potentials, parameterized by w 
and n, signed as (a) — (f) in Fig. |H]and drawn in Fig. have been shown in Fig. |5J 



n 




FIG. 3: Domains of parameters on the plane (w,n) for which the Cardassian models described by 
the system of differential equations (jl9|) start to accelerate for a > 0. In these cases the potential 
function V{a) of the system has a maximum for a > and is upper convex. 
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FIG. 4: Plot of the potential functions describing Cardassian evolution of the universe for different 
parameters (w,n) chosen and signed in Fig. |3]as (a)-(f). 
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FIG. 5: The phase portraits for the Cardassian potentials drawn in Fig. 01 All representative 
phase portraits for different choice of pair (w, n) are shown. The observation of SNIa favors the 
cases (a),(b),(c) for which the static critical points are present. In the neighborhood of this 
solution there are trajectories with very slow expansion or loitering [lj]]. In the case (d), there 
is a complex critical point at (0,0) and for the cases (e),(f) there is no critical points at a finite 
domain. The phase portraits (a), (b) and (c) are topologically equivalent. The critical point at 
infinity (x = oo,x = oo) is a stable node for expanding models. It represents a big-rip singularity in 
the future while x = 0, x = represents the initial singularity which is dominated by the standard 
matter term. Note that models with a big-rip singularity in the future are typical because expanding 
models have a global attractor at the circle at infinity. 

Main advantage of this reduction procedure is the application of dynamical system ap- 
proach to the case of multi-fluids and in general to any fluid for which energy density can 
be expressed as a function of scale factor. For illustration let us consider a universe filled 
with matter and radiation. Then p = p mj o a_S + Pr,o a_4 and an additional term named the 
Cardassian term can be interpreted as some fictitious fluid with energy density pde for which 
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the equation of state is pde = wdePbe, where 

1 



u>de = (n - 1) + 



1 + aVt m j)/Vt r Q 

In this case p c g = p + pde and the potential function for the corresponding Hamiltonian 
system takes the exact form (see Table HJ). 

Let us now consider a class of the so-called generalized Cardassian models, given in the 
form 



x 2 



~2=9(^ (25) 

where x is a dimensionless scale factor expressed in terms of the present value of the scale 
factor ao; x = a/ao; dot denotes here differentiation with respect to the new time variable r\ 
such that dt/Ho = drj. 

Note that function g(x) plays the role of effective energy density in standard Cardassian 
cosmology. We rewrite the generalized Friedmann equation to the form 

£ = Y-9(x)j = 0. (26) 

Equation is an energy function, while —g(x)x 2 /2 plays the role of potential energy for 
particle-universe of unit mass moving in one-dimensional potential. 

Differentiation of (j2*U|) with respect to t] gives the Newtonian equation of motion in the 
form 

dV x 2 
x = --7^7 and V(x) = -g(x) — . (27) 

The Hamiltonian for our problem takes the following form 

H= l -P x 2 + V(x) ) (28) 

where P x = x. 

Equation (|27j) can be reduced to the form of system of second order differential equation 

x = y, 

with first integral x 2 /2 + V(x) = 0. 

It is interesting that all knowledge that is needed about the stationary solution (critical 
points) and its stability can be simply deduced from the geometry of the potential function. 
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TABLE I: The Cardassian models 



The model 


Form of the potential function 


References 


Cardassian model with 
both dust and radiation 




[20] 


Generalized Chaplygin 
gas as a Cardassian model 


V(x) = -I 


-1 , r, 2^ „ , 1"A 3 V /(1 + Q> 
S2 m ,oa: +S2chiip,0^ 1 A » + a 3(l+«) ) 


121. 221 


Cosmology arises from 
MOND gravity 




fim.O^" 1 + n MOND,0 logo! + (1 — £2 m ,o) 


123. 241 



(a) (b) 




0.2 0.4 0.6 0.8 1 1.2 1.4 0.2 0.4 0.6 0.8 1 1.2 1.4 




FIG. 6: The phase portraits for the Cardassian potential functions from Table [I] (a), (b), (c)- 
Cardassian models with dust and radiation for n = —1,0,1 respectively, (d)-MOND model. The 
shaded region is the region of accelerated expansion of the universe. The systems (a), (b), (c) 
are structurally stable because of the presence of saddle points in a finite region. 

We can find how different solutions are distributed on the phase plane. Whether they 
are generic (typical) or exceptional (of zero measure). In terms of dynamical system notions 
the question is whether they constitute a global attractor or a separatrix. It is interesting to 
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check, for example, under which conditions as an asymptotic state for late time we obtain de 
Sitter solution-exit for inflation (acceleration). We can prove that, if the potential function 
for large x (x — > oo) is negative and it is a homogeneous function (of its argument) of degree 
two, then as a critical point we obtain stable critical point at infinity. To prove this it is 
sufficiently to consider the system in projective coordinates z = 1/x, u = y/x. Hence we 
obtain the interesting critical point z = 0, u = uo (it represents de Sitter solution) if only 
V(z) satisfy the condition dV/dz = 2/zV(z). Therefore V(z) oc z 2 . From the linearization 
matrix and the characteristic equation we find the eigenvalues of the linearization matrix 
Ai = —uq, \ 2 = —2uq. They are negative for expanding models. Because both eigenvalues 
are real and negative, the critical point is a stable node. It is located on the intersection of a 
circle at infinity z = — oo < u < +00 and the trajectory of the flat model u 2 /2z 2 +V(z) = 0. 

The potential function for different generalizations of Cardassian models are shown in 
Table Q| 

III. EXACT SOLUTION AND DUALITY RELATION 

For our investigations of exact solutions and transformation relations between FRW model 
solutions with the sub-negative equation of state (1 + 7 > 0, p = jp, 7 = const.) and such 
models with super-negative equation of state (7 < —1), it can be useful to rewrite the basic 
dynamical equation to the form 



where x = a/a is the scale factor in the units of its present value a ; S\ = pi^/iZHo) and 
^Card.o = 1 — ^t,o are the density parameters for matter which satisfy the equation of state 
p = jip (7 = const.) and for some fictitious fluid which mimic the Cardassian term of p n 
type; Qk,o — for simplicity. 

The equation of state of the fictitious fluid which reproduces the Cardassian term is 



Therefore, if n(l + 7) < 0, then an additional noninteracting fluid possesses property of 
super-negative pressure (i.e. phantom matter). The generic property of phantom cosmology 
is the existence of a new type singularities, since they appear for infinite values of the 




(30) 




(31) 



1(3 



scale factor. Hence, if we consider a standard matter (which doesn't violate weak energy 
condition) then it is sufficiently that n < for the appearance of big-rip type singularity. 
In this case the Cardassian term (fluid) dominates for a large x (a ^> a ) while the material 
term is important near the initial singularity. If the Cardassian term dominates the material 
term we obtain the following asymptotic 

x(t) oc (t - t )*&+V. (32) 

From (1321) it is clear that the big-rip singularity appears in a finite time to for Cardassian 
models if only (a) n < and (1 + 7) > or (b) n > and (1 + 7) < 0. In the last case 
x(t) oc (t— 1 ) 2 ^ 3< - 1+7 ^ for large x and Cardassian term is negligible (pure phantom cosmology 
approximation is obvious). 

To find the exact solutions of (J30|) it is useful to introduce a new time variable r\ instead 
of original one t, namely 

3(1+7) 

t — > rj: x 2 dt = drj. (33) 

Of course the new time variable is a monotonic function of the original time and exact 
solutions can be expressed in terms of rj in the following way 

-2/3 



x{rj) 



(34) 



where (3= [3(1 + 7)(1 + n)] _1 . 

In the special cases of domination of the Cardassian contribution and matter we obtain 
the following asymptotic solutions near the big-rip singularities 

^i,0 = 0) ^Card,0 = 1) 

2 

x{jf) oc exp [rj — i]o)t x(t) oc (t — t ) ^"(h-t) (35) 

and 

^,0 = 1) ^Card,0 = 0) 

2 

x(rj) oc exp (rj — rj ) , x(t) oc (t — t ) 3(1+7) • (36) 

Let us note that in time i] the big-rip singularity is moved to infinity. Finally, we can 
conclude that analogical type of unwanted future singularity when infinite energy density is 
reached during a finite cosmological time appear in Cardassian scenario. Moreover it can 
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appear without assumption of super-negative pressure for matter. The existence of this kind 
of singularity is a generic property of the Cardassian models when n < 0. 

It is easy to notice that equation (J3U|) . as well as relation (J33|) . preserves its form under 
the change both position variables and sign of equation of state parameter (1 + 7) 

and (I+7) -> -(1 + 7)- (37) 

From this observation we obtain a simple rule which gives us a solution of the Cardassian 
models in phantom's epochs if we know a solution in quintessence epoch with sub-negative 
equation of state. If x 7 is the solution of (f3"UJ) then xZ^ 1+2 ) * s a ^ so the solution of (jHSJ)- 
In other words, taking 7 = —2/3 like for topological defects we immediately obtain as a 
solution x~ l for phantoms (7 = —4/3). 

Let us now briefly comment this property. The existence of symmetry relation (|37|) 
between the phantom models and quintessence one indicates the presence of duality in the 
dynamical behavior of the phantom and the quintessence field alike to the scale factor duality 
symmetry in super-string theory [25I l^. 

IV. RECONSTRUCTION OF THE DYNAMICS FOR THE CARDASSIAN 
MODEL FROM THE SAMPLE OF KNOP'S SUPERNOVAE DATA 



It was shown in our previous papers |7j, that the potential function V[a;(^)], which 
completely determines the dynamics of the Cardassian model (see section II B), can be 
simply reconstructed from the observed magnitude-redshift relation of distant supernovae. 

As it is well known the relation between the apparent magnitude m and the dimensionless 
luminosity distance T>l is 

m = M + 51og[P L (^m,o,^Card,o,^, z)], (38) 

where M. = M — 5logH + 25; M is the absolute magnitude; H is the Hubble constant 
and the dimensionless luminosity distance for the flat model takes the form 

T^L^mfi, ^Card,0, n, z) = H dL(Q m fl, ^Card.O, H , z) 

, f dz' M , f z dz' 

= #0 1 + z) / -— — - = (1 + z 



o H(n mfi ,n c ^n,H 0j z') J y/n mfi (i + V) + n C ard,o(i + z') 3n ~ 2 ' 

(39) 
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Due to the existence of such a relation it is possible to calculate the potential function V 
which is 

1 



V[a(z), fi m)0) ^card.o, n, H , z\ = -^H 2 (Vt mfl , ^card.o, n, H , z)a 2 (z) 



d ( dL(Q miQ ,QcaTd,o,n, H , z) 
dz 



1 + 2 



a 2 (z) 



—Hi 
2 



^m,o(l + z ) + ^Card,d(l + z 



,3n-2 



(40) 



-3n+2 



(41) 



Relation f|4(J|) can be rewritten in the following (dimensionless) form 

V[x(z), n m>0 , fi C ard,0, n] = -- QmfiX' 1 + ^Card 

We determine the model parameters (S- m ,o, i- Card , n, M) as well as the potential function 
V using a x 2 minimization procedure, where x 2 is described by the formula 



X — a 2 ' v 4 ^v 

where m j is the measured magnitude; m t j can be calculated from the equations (|38j) 
and for given values of parameters f2 m> o, ^Card.o, n, M- and z\ a m j is the magnitude 
measurement uncertainty. The summation in the equation ()42|) is over all of the observed 
super novae. 



TABLE II: Results of the statistical analysis for the Cardassian model from distant type la super- 
novae data (Subset 1 of 58 supernovae compiled by Knop et al. [271]). obtained from the best fit 
method (x 2 minimization) and from the likelihood method (denoted with max(£)). In the case (a) 
the analysis was prepared with fixed A4 (chosen as the best-fitted value of Ai for the model). The 
probability density distributions (from the likelihood procedure) for the estimated cosmological 
parameters are shown in Fig. |HJ 





^m,0 


^Card,0 


n 


M 


x 2 


Zmax 


Vmax 


best fit 
max(£) 


0.55 

n r; 7 +0.05 


0.45 

n 4S+ - 05 

u - 4J -0.05 


-4.46 
-4.36±i;f? 


-3.62 
-3.62 


53.42 


0.189 


-1492.96 


best fit 
max(£) 


0.55 
^4+ - 07 

u - O4 -0.06 


0.45 
46+ 06 


-4.46 

_4 1 4+3.55 


-3.62 

_o fil +0.08 


53.42 


0.189 


-1492.96 



Assuming that the measurement uncertainties are Gaussian we can simply calculate the 
likelihood function £ from a chi-squared statistics in the following way 

C = 1 exp ( - ^ j , N = JJJJ exp ( - |- j dn mfi dQ c ^,odMdn, (43) 
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FIG. 7: Confidence levels on the plane (£l m ,o, n) for the Cardassian model obtained by the integra- 
tion of the likelihood function £(Q mj o, f2card,o> n ) over ( a ) ^m,o> -M. was chosen as the best-fitted 
value of M for the model in minimization procedure; (b) f2 m .o an d M for the subset 1 of 58 
supernovae compiled by Knop et al. j22j. -M. was chosen as the best-fitted (see Table Hlj) . 



where N is the normalization coefficient given as an integrate of the likelihood function over 
all probability variables. 

We use the "full primary subset" (Subset 1) of 58 supernovae compiled by Knop et al. 
3| for the statistical analysis. Both the stretch correction and the host-galaxy extinction 
correction were applied to the sample At first the \ 2 minimization procedure was 

prepared for the sample. Results of the analysis are presented in the Table |H] Then the 
three-dimensional space of estimated parameters (fi m ,o, ^Card,o, n) was divided into a grid 
and at each point of grid the likelihood function was calculated for previously best-fitted 
value of M. parameter. It is important to cover the entire space of estimated parameters 
where the likelihood function differs from zero. In our analysis the exploded space included 
ranges of parameters: Q m ,o G [0 : 1], f2c ar d,o G [0 : 1], n G [—13 : 1]. 

The confidence levels for pairs (fi m ,o, n), obtained through the integration of the likelihood 
function £(fi mi0 , ^card,0) n) over the fic a rd,o variable, are drawn in Fig. [3(a). 

The probability density distributions for cosmo logical parameters Aj (fi m0 , fi Card0 , 
n) from the likelihood analysis, obtained by the integration of the likelihood function 



20 




FIG. 8: The probability density distributions for cosmological parameters Ai (£l m ,o, ^Card,o> 
n, M) from the likelihood analysis obtained by the integration of the likelihood function 
£(^m,o> ^Card,o> n i -M) over remaining probability variables C{Ai) = Yik=£i dAkC(Aj). The 
cases (a), (c) and (e) correspond to the analysis with fixed A4 = —3.62 
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FIG. 9: Confidence levels on the plane (f2 m> o,n) for the flat Cardassian model and different values 
of M parameters (a) M = 3.62 (best fitted value for the model), (b) M = 3.58 and (c) M = 3.54 
(best fitted value for curved model with non-relativistic matter and cosmological constant). Note 
that results of the statistical analysis strongly depend on the value of M parameter. 

As it is shown in Fig.|U]the likelihood function strongly depend on the value of A4 param- 
eter so it is very important to choose it as a best-fitted value for the model or (for a more 
robust determination) to treat M. as an additional probability variable and integrate full like- 
lihood function also over M.. The analysis with marginalization over M. was prepared in the 
similar way. The results of the analysis were presented in the Table [0] The confidence levels 
on the plane (Q m Q ,n) and one-dimensional probability density distributions for the cosmo- 
logical parameters fl m ,o, ^card,o ; n > -M are shown in Fig. [3(b) and Fig. |HKb),(d),(f),(g) 
respectively. Of course, in this case, the uncertainties of the cosmological parameters are 
bigger because of one additional estimated parameter. 

It is well known that there is a simple interpretation of the Cardassian term as an ad- 
ditional noninteracting fluid with pressure p = [n(l + w) — l]p, where the standard term 
is obtained from the fluid p = wp and w = const. In the special case of dust we obtain 
p — (n — l)p. Let us note (see Fig. EK a )) that the supernovae data favor the high density 
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matter model of the universe with n < at the 99% confidence level (it is not truth for the 
cases (b) and (c)). From the previously presented interpretation of the Cardassian model 
as a FRW model filled with a mixture of noninteracting fluids we can conclude that the 
universe is now dominated by the mater with a super-negative p_ressure p < —p. At first the 
fact was noticed by T. Roy Choudhury and T. Padmanabhan 



(a) 



(b) 




V-reconstruction from 58 supernovae 
compiled by Knop at al. (2003) 




V-reconstruction from 58 supernovae 
compiled by Knop at al. (2003) 



FIG. 10: Confidence levels (la and 2a) for the potential function V{z) given by the formula (|41|) 
obtained from the likelihood method (see Table Hlj) . The case (a) corresponds to the analysis with 
fixed M = —3.62. The best fit is represented by the solid black line. 




FIG. 11: Phase space for the system (|19f) . where the potential function F[x(^)] given by the 
equation (fFTj) has been reconstructed from the SNIa observational data [22] ■ The shaded domain 
of phase space is the domain of accelerated expansion of the universe. The dashed line represents 
the flat model trajectory which separates the regions with negative and positive curvature. 



We learn from the previous dynamical systems analysis that all what is required for unique 
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determination of the Cardassian models contain the diagram of the potential function. Since 
the Hubble function is related to the luminosity distance (for the flat model) by the simple 
relation 

it is possible to determine both the potential function for the system (or equivalent p e g(z)) 
and then quintessence parameter w(z) = p c s/p c a- In our previous papers [?], Q it was 
demonstrated that the reconstruction of the potential function can be performed for the 
general class of FRW models with an equation of state in the form p = w(a)p. Such 
an approach to the reconstruction of the dynamics, in terms of the potential (or energy 
density), directly from observations in the model independent manner has two advantages: 

1) the dynamics can be derived only from the single function and the description has a 
natural interpretation - particle-like description; 

2) the reconstruction of V[a(z)] is a simpler procedure then the probe of the dark energy 
in term of quintessence parameter w(z), for reasons pointed out by Maor, Brunstein and 
Steinhardt j2jj|. 

Because of many problems in determining the quintessence parameter w(z) we propose 
the reconstruction of V[a(z)] simpler alternative to description the dynamics 

of quintessential cosmology. Later some authors [soj], for the same reasons, use a rescalled 
quantity of the potential, i.e., p(z) = —6V(z)(l + z) 2 having the sense of effective energy 
density. Advantage of use of the potential function is that the dynamics can be completely 
determined. Our paper shows the advantage of using V(z) instead of w(z) in investigation 
of the dynamics of quintessential models. 

Note that both deceleration parameter q and quintessential parameter w are simply re- 
lated to the elasticity of the potential with respect to the scale factor iy(a). 

q = --I v {a), w(t) = --(I v (a) + l), 

w = --l'v(a), (45) 

where Iy = dlnV/dlna is the elasticity parameter. 

Because trajectories of the system lie on the zero energy level (Hamiltonian constraint) 
we obtain 

d I d L (z) 



V(z) = ~(l + z) 



dz \ 1 




(46) 
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where di{z) is fitted by using the Cardassian model itself rather then the polynomial ansatz 
(as in our previous work) because it gives better results. 

Therefore from the above formula we can simply obtain the effective energy density 



Note that in both last formulas only first derivatives of di(z) appear while w(z) depends 
on the cIl(z) through the second order differentiation. 

In the Fig. ^2 we present the best fit for the reconstructed potential function of the 
Cardassian scenario and the confidence levels (la and 2a) for the potential obtained from 
the likelihood statistical analysis while the phase portrait for the reconstructed potential 
function is shown in Fig. El In the statistical analysis we consider the polynomial function 
and the power low function. The lowest x 2 was obtained for the latter so we assume the 
power low ansatz for the reconstruction of the potential function. 

It is very interesting that from the reconstructed potential we can immediately obtain the 
phase portrait. It is a consequence of the fact that the potential uniquely determines the dy- 
namics of the Hamiltonian flow. This portrait is representing the evolutional path of the Car- 
dassian model for all initial conditions. Note that the phase portraits from Fig. |HKa),(b),(c) 
are topologically equivalent to that obtained from the reconstruction (see Fig. fTTj) . 

V. ON SOME OTHER SPECIAL FEATURES OF CARDASSIAN MODELS 

It seems to be reasonable to require the phase portrait equivalence to be fulfilled. It means 
that any model of the universe should describe its dynamics without contradiction with 
astrophysical data. Therefore, the acceptable model should have phase portrait (qualitative 
dynamics) equivalent to the phase portrait reconstructed from observations. 

In this section we would like to present how Cardassian models can explain some cosmo- 
logical problems j^. The Cardassian models are able to provide a natural explanation for 
our observations of the universe which is almost isotropic and flat. In standard FRW cos- 
mology with perfect fluid we can observe that at large scale factor a the curvature term falls 
off faster than a -3 ( 1+1u ) term whenever w < —1/3, i.e., some special kind of fluid is required 
for explanation of the flatness problem. By contrast, in the Cardassian scenario the flatness 
problem is naturally solved because the ka~ 2 term falls off faster than the Ba~ 3( - 1+W ^ n term 




(47) 
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provided n < 2/(3(1 + w)). Therefore the period of dust-dominated evolution of the uni- 
verse with n < 2/3, which satisfies following Cardassian assumption, can solve the classical 
flatness problem. 

Let us consider the class of anisotropic but homogeneous cosmological models, namely 
Bianchi model with perfect fluid, at late time of evolution. In these models anisotropy 
rate arises from scalar anisotropy, which falls off no slower than a 2 oc a~ 2 at late time, 
where a is an average scale factor. In the homogeneous, but anisotropic flat Cardassian 
cosmological model the basic equation is modified by addition of a a 1 term to its right-hand 
side. Therefore, the isotropy problem is solved by accelerating universe with < w < 2/3 
because scalar anisotropy a 2 falls off faster than the a~ 3 ^ 1+w ^ term. If we consider Cardassian 
model with n < 2/(3(1 + w)) then the anisotropy term can not dominate the expansion. 
Thus the anisotropy problem can be solved in a dust-dominated universe by an interval of 
Cardassian evolution if n < 2/3. Of course, in the standard cosmology if p + 3p < (strong 
energy condition is violated) the matter term falls off more slowly than a 2 as t — > oo, i.e., 
anisotropy is damping during the accelerating phase of the evolution whereas the Cardassian 
scenario offers a natural explanation of anisotropy problem without assumption of a special 
form of dark matter. It is a consequence that Cardassian models solve the flatness problem. 
Therefore these models do not offer a new solution of the isotropy problem. This situation 
in some sense is analogical to that which take place in the context of VSL cosmologies 3^| • 
However let us note that VSL models in contrast to the Cardassian models can solve the 
horizon problem in the past. 

Let us now consider some general situation in which we have w = w[a(t)]. Such a case 
appear if we consider both matter and radiation for example. Then the following theorem 
establish the sufficiently condition for the presence of particle horizon. 

When all events whose coordinates in the past are located beyond some distance dn then 
can never communicate with the observer at the coordinate r = (in R-W metric). We can 
define the distance du as the past event horizon distance 

d H (t) = a(t) / — -c = a(t)I. (48) 

Of course, whenever / diverges as to ~* 0, there is no past event horizon in the spacetime 
geometry. Then it is in principle possible to receive signals from sufficiently early universe 
from any comoving particle like a typical galaxy. If the if integral converges for t — > then 
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our communication with observer at r = is limited by what Rindler has called a particle 
horizon. The particle horizon will be present if energy density is growing faster then a~ 2 ~ e 
as a — > (e > 0). Therefore if d > a~ e / 2 (d > A) then there is a particle horizon in the past. 

Our investigation of the particle horizon here is independent of any specific assumption 
about the behavior of a(t) near the singularity or of specific form of the equation of state. 
If one assumes a linear equation of state p = wp and w = const., then Friedman's equations 

2 

imply the following behavior for a(t) ±2 (£) 3 (t +1 ) near the singularity t — 0. The integral / 
would thus diverge only if 7 < —1/3, i.e., only if the pressure p becomes negative. This is 
the condition for solving the horizon problem and it is identical to that for the solution of 
the flatness problem (p/3 term will dominate the curvature term in a long time evolution). 
We can also show here that the integral J would diverge only if the pressure of the cos- 
mological fluid takes negative values in the general case of w(a(t)). Conservation condition 
can be rewritten in the form 

a 3 dp = d[pa 3 (l + w(a))}. (49) 

We can verify that boundedness of a(t) means also that pa 2 remains bounded near the 
singularity. Therefore, pa 3 — > as t — > 0. By integrating both sides of the equation written 
above from to t we obtain —3 j^w{a)pa 2 adt = a 3 p(t) ^ 0. 

Consequently, w(a) must assume a negative value without any specific assumption about 
the equation of state. 

VI. CONCLUSIONS 

In this paper we apply the qualitative cosmology analysis to the general class of Car- 
dassian models which have become rather an alternative to the dark energy models. We 
show that while the application of qualitative methods allows to reveal some unexpected 
properties like the big-rip singularities, the SNIa data uniquely reconstruct the phase plane 
in the model independent way. The reconstruction of this phase plane should be treated as 
a necessary condition for any model which wants to explain the present acceleration of the 
universe. 

Theory of dynamical systems which offer a possibility of investigating the space of all 
solutions for all admissible initial conditions is used in analysis of Cardassian models. We 
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demonstrate a simple method of reducing of the dynamics to certain two-dimensional dy- 
namical system. One of the features of this reduction is the possibility of representing the 
model as a Hamiltonian system in which the properties of the potential function can serve 
as a tool for qualitative classification of all possible evolution scenarios. It is shown that 
some important features like resolution of the flatness problem or horizon problem and its 
degree of generality can be visualized as domains on the phase plane. Then one is able to see 
how large is the class of solutions (labelled by the initial conditions) leading to the desired 
property, i.e., this class is generic or non- typical. 

Applying the techniques of the particle-like description, developed by us earlier, to a new 
data set from the Knop's sample we show that the Cardassian model with n < is favored 
by the data on the 99% confidence level. On the other hand it is just the case of appearance 
of unexpected big-rip singularities in the future evolution of the model. 

Although the developed formalism is mainly adopted to the investigation of the Cardas- 
sian models it can be useful to include additionally dissipative effects of bulk viscosity to 
the model. Moreover for these class of models de Sitter solution is admissible as a global 
attractor in the future. 

The main results are the following: 

1. We show the effectiveness of using of dynamical system methods in analysis of the 
wide class of Cardassian models and their generalizations. 

2. The unwanted big-rip singularities at a finite time were detected as a generic for the 
case n < (a global attractor in the future). 

3. The exact solutions and some duality relations were found. 

4. Additionally we showed that the existence of this kind of singularities is acceptable on 
the 99% confidence level from Knop's sample of SNIa data. 

5. The particle-like approach was adopted to describe the dynamics as well as the po- 
tential function was reconstructed from the observations. In the fitting procedure the 
power law function are used because it gives the lowest value of x 2 . 

6. Because of interpretation of the Cardassian term as an effect of additional fluid one 
can claim that if FRW equation is "correct" then phantom fields are required at the 
99% confidence level. 
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Finally, our study showed that the Cardassian model with n < is strongly favored by 
the Knop's SNIa data (without any priors on matter content). But, this model predicts 
rather unexpected future of the universe - a big-rip singularity. However, the estimation of 
characteristic time of the big-rip singularity will be a subject of our next paper. 
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